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RANDOM  VARIANCE  OPTION  PRICING: 
EMPIRICAL  TESTS  OF  THE  MODEL  AND  DELTA-SIGMA  HEDGING 

Abstract 


We  present  an  empirical  analysis  of  the  random  variance  option 
pricing  model  by  studying  prices  on  call  options  from  the  CBOE.   The 
random  variance  model  is  compared  with  the  Black-Scholes  model  and 
both  models  produce  prices  that  fit  actual  option  prices  quite  well. 
For  both  models  we  find  pricing  biases  with  respect  to  the  strike 
price,  time-to-maturity,  and  volatility.   There  are  consistent  strike 
price  biases  and  tirae-to-maturity  biases  in  our  sample  for  the  Black- 
Scholes  model.   Similar  biases  are  present  in  the  random  variance 
model,  but  the  direction  is  mixed  (both  underpricing  and  overpricing). 
We  examine  an  alternative  hedging  strategy  with  options  in  which  we 
hedge  both  the  price  risk  and  the  volatility  risk.   We  find  that  over 
two-day  trading  intervals  one  can  reduce  the  variability  of  a  hedged 
option  position  by  hedging  both  price  and  volatility  risk. 


RANDOM  VARIANCE  OPTION  PRICING: 
EMPIRICAL  TESTS  OF  THE  MODEL  AND  DELTA-SIGMA  HEDGING 


I.   Introduction 

Recent  papers  by  Hull  and  White  (1987a),  Scott  (1987),  and  Wiggins 
(1987)  have  examined  the  pricing  of  options  on  assets  with  stochastic 
volatility.   The  general  approach  to  pricing  options  in  these  papers 
is  to  apply  the  methods  of  Cox,  Ingersoll,  and  Ross  (1985a)  and  treat 
volatility  as  a  random  state  variable.   The  papers  by  Scott  and  Wiggins 
contain  some  preliminary  empirical  analysis  of  these  models.   In  this 
paper,  we  present  a  more  complete  empirical  study  of  the  random  vari- 
ance model  by  studying  transaction  prices  from  the  CBOE.   In  the 
second  section  of  the  paper,  we  present  a  method  for  estimating  the 
parameters  of  the  volatility  process  for  the  underlying  stock.   These 
estimates  are  needed  in  order  to  calculate  option  prices  with  the 
model.   In  the  third  section,  we  examine  various  option  pricing 
issues.   We  price  options  taken  from  the  Berkeley  Options  database  for 
the  year  1983  and  we  compare  the  performance  of  the  random  variance 
model  with  the  Black-Scholes  model.   Our  analysis  includes  an  examina- 
tion of  three  option  pricing  biases  that  have  been  documented  in  the 
literature:   the  strike  price  bias,  the  time  to  maturity  bias,  and  the 
volatility  bias.   In  Section  IV,  we  examine  the  hedging  implications 
of  the  random  variance  model.   In  the  random  variance  model,  one  must 
hedge  both  price  and  volatility  risk  in  order  to  eliminate  the  risk  of 
an  option.   To  test  the  importance  of  volatility  risk,  we  use  the 
models  and  actual  prices  to  form  hedged  positions,  and  then  we  calcu- 
late the  variability  of  these  hedged  positions.   To  hedge  both  price 
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and  volatility  risk.,  one  must  hedge  an  option  with  the  stock  and 
another  option  or  hedge  an  option  with  two  other  options.   We  follow 
Hull  and  White  (1987b)  and  refer  to  this  hedging  strategy  as  delta- 
sigma  hedging.   We  compare  the  variability  of  the  delta-sigma  hedge 
with  the  variability  of  the  familiar  delta  hedge,  which  is  used  by 
option  traders. 

We  begin  by  reviewing  the  theory  of  random  variance  option  pricing 
applied  to  call  options  on  stocks.   We  use  a  diffusion  model  that  is 
similar  to  the  models  examined  in  Hull  and  White,  Scott,  and  Wiggins: 

dP 

—  =  adt  +  adz 

(1) 

—  =  f(a)dt  +  ydz. 
a  L 

where  E(dz  ,dz~)  =  5dt  and  the  interest  rate  is  a  constant  r.   We  con- 
sider specifically  an  Ornstein-Uhlenbeck  (mean  reverting)  process  for 

—  12     — 

Ina   so  that  d2.no  =  8(a-2.na)dt  +  ydz2  and  f(o)  =  -ry     +    S(a-2na).   The 

fundamental  partial  differential  equation  (P.D.E.)  for  the  price  of  a 
call  option  written  on  the  asset  with  price  P  is 

-H3  +  \  HU02P2  +  Hl25Ya2P  +  \  H^yV  -  Hr 

+  H^P  +  H2[y  y2  +  B(I-£na)  -  X*]  =  0, 

where  X*  is  the  risk  premium  associated  with  volatility. 

The  solution  to  the  P.D.E.  and  the  boundary  conditions  for  call 
options  has  the  following  form: 

H(Pt,at,t,T)  =  E(e"r(T"t)max(0,PT-X}jPt,at),  (2) 
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where  X  is  the  strike  price,  t  is  current  time,  T  is  the  expiration 
time,  and  E  indicates  the  expectation  taken  with  respect  to  the  risk 
adjusted  processes: 

dP 

■5—  =  rdt  +  adz 


(3) 


—  =  (f(a)-X*)dt  +  ydz0. 


Hull  and  White  and  Scott  have  shown  that  if  stock  price  changes  and 
volatility  changes  are  uncorrelated  (5=0),  then  the  solution  can  be 
simplified  to 

CO 

H(P  .a  .t,T)  =/   [P  N(d1)-Xe"r(T"t)N(d_)]dF(V),       (4) 
t   t         ~    t    i  I 

T 

t  2 

where  V  =  f  a,  xds  and  F  is  the  distribution  function  of  V.   N(«)  is 

t  (s) 

the  cumulative  standard  normal  distribution  function  and  d  and  d„  are 

the  same  d..  and  d  used  in  the  Black -Scholes  model,  except  V  is  used 

2 
in  place  of  a  (T-t).   The  solution  in  equation  (4)  is  less  expensive 

and  easier  to  use  in  an  empirical  study  and  we  find  in  the  next  sec- 

2 
tion  that  estimates  of  6  have  very  large  standard  errors.    We  calcu- 
late prices  by  Monte  Carlo  simulation  of  equation  (4),  and  we  use  the 
antithetical  variate  method,  which  is  described  in  Boyle  (1977). 


II.   Parameter  Estimates  of  the  Volatility  Processes 

To  implement  the  random  variance  model  in  either  equation  (2)  or 
equation  (4),  one  must  have  values  for  the  risk  premium  and  the  param- 
eters of  the  volatility  process.   In  this  section  we  discuss  a  method 
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of  moments  estimator  for  the  parameters  of  the  volatility  process. 
Estimating  a  value  for  the  risk  premium  is  discussed  in  the  next  sec- 
tion.  Our  approach  is  to  apply  a  discrete  approximation  of  the  model 
to  daily  stock  returns  and  then  estimate  the  parameter  values  by 
applying  the  method  of  moments  as  discussed  in  Hansen  (1982).   We  use 
one-day  intervals  for  both  the  parameter  estimation  and  the  Monte 
Carlo  simulation  of  the  model.   We  start  with  the  price  solution  for 
the  model  in  (1): 


t      ,   «       t 

P   =  ?  exp{/   (ct  -  j   a^)ds  +  /   a  dz  (s)}. 
tu      Q      ^  s       0s 

We  use  the  following  approximation  for  AilnP   over  a  short  time 

,  3 
interval : 


A*nPt  =  (a  "J  Vl)At  +  at-lAZ!t' 


which  we  rewrite  as 


A£nPt  =  uAt  +  0-  ^A*   . 

1   2 
where  u  =  a  -  —  a  _1 •   Since  Ina   is  an  Ornstein-Uhlenbeck  process,  we 

have  the  following  equation  for  ina     at  discrete  points  in  time: 

Una     ■  e   Jtna  ,  +  a(l-e   )  +  e 
t         t-1  t 

where  e   is  normally  distributed  with  mean  zero  and  variance 
t 

2    -26 
Y  (1-e  '  )/28.   Letting  At  =  1,  we  have  the  following  model  for  our 

empirical  analysis: 
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A*nPt  =  V  +  VlUt 

2,no\  =  a  +  pina    +   e 
t  t-1    t 

where  u  is  a  standard  normal  random  variable  and  e  is  normally  dis- 

2 
tributed  with  mean  zero  and  variance  a   .   We  allow  correlation  between 

e 

stock  price  changes  and  volatility,  E(u  e  )  ■  6a  ,  and  we  assume  that 

E(u  u  ,  )  =  0  and  E(e  e  ,  )  ■  0  for  k  *■   0.   The  parameters  ~a,    6,  and  y 
t  t-k  t  t-k 

2 
can  be  easily  recovered  from  the  estimates  of  a,  p,  and  a  as  follows: 

—    a 
a  = 


1-P 


=  -2,np 


2    2  (-2 top) 

Y  -  a  — 

£  (1-P  ) 


Now  consider  the  estimation  of  the  fixed  parameters,  a,  p,  a  ,  and 
6.   Maximum  likelihood  estimation  of  these  parameters  is  intractable 
because  we  cannot  derive  a  closed  form  for  the  distribution  of  a 
sample  of  time  series  observations  on  stock  returns.   For  this  reason, 
we  employ  a  method  of  moments  estimator.   Our  approach  is  to  use  a 
variety  of  moments  to  identify  and  estimate  the  four  fixed  parameters. 
Since  we  are  not  interested  in  estimating  u,  we  work  with  deviations 
from  a  standard  estimator  such  as  the  sample  mean: 


xt  E  A*nPt  -  |i  =  a^iy 


The  second  and  fourth  moments  of  x  are 


E(xj:)  =  exp{2(-4-)  +  2(— £=■)} 

1-P 
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2 

c 

)  +  2C— 
•p 


2 

.  o 

EU   )    -  3  exp{4(-^-)  +  8(— ^) } 
C  p      1-p 


From  sample  estimates  of  the  two  moments  we  can  compute  estimates  for 

a2 
a  v    .  „  e  „    ,  a 

■P- 
that  kurtosis  for  the  stock  return  is  a  function  of  ( ~) : 


(- )  and  ( tt) .      (- )  is  the  mean  reverting  level  for  Una.      Note 

1-P       !_p2     1-P  ^ 

^P1 


E(x  )  a2 

K  = 2~T  =  3  ex?f4( ?>}» 

[E(xJP  1-P 


and  we  have 

-2 


( 2")  =  "4  in^' 

1-p 


If  the  sample  kurtosis  is  not  greater  than  three,  our  parameter  esti- 

<y- 

p 

mation  breaks  down.   We  are  essentially  setting  ( ~-)  to  explain  the 

estimated  level  of  excess  kurtosis  in  stock  returns.   For  (- ),  we 

1-p 

have 

A* 

A. 

(-; )  -  TT  Jtn(E(x  ))  -  y   infcr). 

1-p    z       t     *»    J 

We  need  additional  moments  to  identify  estimates  of  p  and  5.   To 
derive  these,  we  examine  the  behavior  of  £n|x  j  ■  ina         +    toju  j, 
where  the  last  term  is  the  log  of  the  absolute  value  of  a  standard 
normal  variate.   From  the  autocovariance  of  £n|x  j,  we  have 
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2 
a 

Cov(fcn|x  j  ,*n|x    j  =  p( j) , 

1-p 

and  from  the  corresponding  sample  moment  we  can  infer  an  estimate  of  p. 

Given  p,  we  can  use  our  previous  sample  estimates  to  compute  estimates 

2 
of  a  and  a   .   We  have  found  that  estimating  p  from  this  set  of  sample 

moments  does  not  produce  a  very  precise  estimate.   Conceptually  one 

could  use  the  entire  autocovariance  function,  or  the  autocorrelation 

function,  for  £njx  j  to  compute  an  estimate  of  p: 

2 

a 

VarUn|x    |  )    =  — %£  +  Var(to|u    |) 
1-p 

2 
.       a 

CovUn|xti,£n|xt_k|)    =  p   (— —■) ,   k   *  0 

1-p 

and   Var(£n|u    J)    =  — .        The   time   series    £n|x    j    is    the   combination   of   a 
to  t 

first    order   autoregressive   process   and   a  white   noise   error    terra,    and 
it   has    the   appearance   of   an  ARMA    (1,1)    model.      Observe    that 

(1-pLHnjxJ    =  a  +  ec_1   +    (1-pL)  to  |u    |  , 

re  L  is  the  lag  operator  (Lx   =  x   , )•   On  the  right  hand  side  we 


whe 


have  a  first  order  moving  average  process,  and  the  following  is  an 
observationally  equivalent  model: 


(l-pLHn|x  j  »  b  +  (l-9L)n 


*«K\  ' b* +  ri^>  \' 
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where  n   is  a  white  noise  error  term.   9  and  Var(n  )  are  functions  of 

a  ,  Var(£n|u  I),  and  p.   The  p  parameter  can  be  estimated  from  the  ARMA 
e        '  t 

(1,1)  model  for  £n|x  j.   We  are  effectively  using  the  entire  autocorre- 
lation function  to  estimate  p  and  the  result  is  a  more  efficient  estimate. 

The  6  parameter,  which  measures  the  correlation  between  stock  price 
changes  and  volatility,  can  be  estimated  from  a  variety  of  moments.   To 

identify  this  parameter  we  must  use  the  cross-covariance  between  x  and 

2  ii 

future  values  of  x  or  future  values  of  £n|x  |.   The  cross-covariance 

function  for  £n|x   and  x  has  the  following  form: 

1    t '  t  ° 

2 

•-1  1      a 

CovUn|x    |,x        )    =  pJ      6a£   exp{(^-)    +  y( — %£ >  } 

J  1-p 

for  j  =  1,  2,  3,  ... 

=  0         for  j  ^  0. 

Jenkins  and  Watts  (1968,  pp.  336-40)  have  noted  that  one  can  improve 

estimates  of  the  cross-covariance  function  by  filtering  the  series. 

If  we  filter  the  series  £n|x  |  by  using  our  estimate  of  p,  we  obtain 

the  following  result  for  the  cross-covariance  function  between  x  and 

yt  =    (l-pL)£n|x  |: 

2 

Cov(yt,x   )  -  $oe   exp{(~-)  +  y(— ^) }  if  j  -  1 

1-p 

=  0  If  j  it  1. 

The  resulting  estimator  is 

5  = expj-^)  -  y  ( ?)J, 

e  1-p 

where  we  use  p  to  compute  y  and  we  compute  the  sample  covariance. 
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To  assess  the  precision  of  our  estimates,  we  calculate  the  variance 
matrices  and  standard  errors  for  the  various  parameter  estimates.   We 
estimate  p  and  9  from  the  ARMA  (1,1)  model  for  £n|x  |  by  using  the 
frequency-domain  method  described  in  Hannan  (1970,  pp.  396-98).   Hannan 
and  Walker  (1964)  both  provide  large-sample  distribution  theories  for 
this  estimator.   The  other  parameter  estimates  are  functions  of  p  and 
the  other  sample  moments.   To  compute  large  sample  variances  for  these 
estimators  we  use  the  following  familiar  result:   for  a  parameter  a 
that  is  a  function  of  other  parameters  or  moments,  f(m),  if  /T(m-m) 
converges  to  a  normal  distribution  with  mean  zero  and  variance  matrix 
V,  then  /T(a-a)  converges  to  a  normal  distribution  with  mean  zero  and 


1  to  Or— )  'V(-x— ).   To  estimate  the  variance  matrix  for  our 
dm  dm 

sample  moments,  we  use  methods  described  in  Hansen  (1982).   Because  our 


variance  equa. 

dm    dm 


time  series  may  be  serially  correlated,  we  use  the  spectral  estimator 
for  this  variance  matrix.    In  addition  to  the  parameter  estimates,  we 
compute  the  sample  kurtosis, 

A   4 
E  (x*) 

K  =  


2   2 

(e  (xpr 


and  its  standard  error  so  that  we  can  test  whether  x   is  normally  dis- 

t 

tributed  (K=3).   In  the  ARMA  (1,1)  model  for  £n|x  |,  we  have  white 
noise  and  an  identification  problem  if  p  =  9.   Hence  a  second  test  is 
based  on  the  t  statistic  for  p  -  9,  and  if  we  were  to  find  that  p  is 
not  significantly  different  from  9  then  we  would  interpret  this  result 
as  implying  p  =  9  =  0. 
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We  use  daily  stock  returns  to  estimate  the  parameters  of  the  vola- 
tility processes,  and  our  data  samples  include  the  NYSE  value  weighted 
portfolio,  the  S&P  500,  and  13  stocks  whose  options  trade  on  the  CBOE. 
We  have  chosen  12  stocks  included  in  the  Dow  Jones  Industrial  Average 
plus  National  Semiconductor  so  that  we  have  one  stock  in  the  sample 
that  does  not  pay  dividends.   The  sample  periods  run  from  November 
1974  to  year-end  1982.   Our  goal  originally  was  to  construct  samples 
of  roughly  ten  years  of  daily  returns  through  1982,  but  we  have  reduced 
the  samples  to  series  with  2048  observations,  which  is  a  multiple  of  2, 
to  take  full  advantage  of  the  fast  Fourier  transform  in  computing  esti- 
mates for  the  ARMA  model.   The  returns  are  computed  cum  dividend  and 
are  taken  from  the  daily  CRSP  tape. 

The  first  step  is  to  compute  x  >  which  is  the  log  of  one  plus  the 
rate  of  return  deviated  from  its  expected  value.   Many  studies  simply 
ignore  the  expected  return  in  daily  return  data  because  the  number 
will  be  quite  small  relative  to  the  daily  price  change  and  the  esti- 
mates are  typically  not  very  precise.   One  method  would  be  to  subtract 
off  the  sample  mean,  but  recent  empirical  research  indicates  that 
there  is  a  strong  day-of-the-week  effect  in  stock  returns.   For  this 

reason  we  have  used  dummy  variables  for  the  days  of  the  week  to  remove 

g 
this  effect  in  calculating  the  x  series  for  each  stock.   Two  of  the 

stocks,  IBM  and  United  Technology,  have  missing  observations.   There 

were  four  missing  observations  and  in  each  case,  we  have  simply  spread 

the  two-day  returns  over  both  days.   In  one  case,  we  were  able  to 

obtain  a  closing  price  on  the  Midwest  Stock  Exchange:   on  January  27, 

1975,  IBM  did  not  trade  on  the  NYSE,  but  after  the  close  in  New  York, 
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trading  resumed  on  the  Midwest  Exchange  and  I3M  closed  up  $20  per 
share. 

The  parameter  estimates  and  test  statistics  are  presented  in  Table 
I.   The  sample  size  for  the  estimator  of  the  ARMA  models  is  2048,  but 
it  is  2047  for  the  other  parameters  and  test  statistics.   The  p  esti- 
mates are  all  close  to  one,  and  most  of  the  estimates  are  within  two 

2 
standard  errors  of  one.   The  resulting  estimates  for  a  and  a     indicate 

e 

plausible  values  and  the  estimated  standard  errors  indicate  a  reason- 
able level  of  precision.   The  standard  errors  for  the  5  estimates  are 
quite  large  and  indicate  unreliable  estimates.   None  of  the  estimates 
are  significantly  different  from  zero  and  the  sample  covariances  from 
which  the  6Ts  are  estimated  are  all  close  to  zero.   The  p  estimates  in 
Table  I  are  consistent  with  estimates  found  in  Poterba  and  Summers 
(1986):   the  mean  half-lives  from  our  estimates  range  from  13  trading 
days  to  a  high  of  268  trading  days.   The  mean  half-lives  for  the  two 
market  indexes,  the  NYSE  and  the  S&P  500,  are  respectively,  70  and  85 
trading  days. 

The  other  statistics  in  Table  I  indicate  rejection  of  the  lognor- 
mal  model  for  stock  returns  and  support  for  the  view  that  there  is 
some  persistence  in  volatility.   The  sample  kurtosis  is  greater  than 
four  for  all  fifteen  series  and  all  of  the  estimates  are  significantly 
greater  than  three,  which  is  the  theoretical  value  if  x  is  normally 
distributed.   The  tests  for  p  =  9  in  the  ARMA  model  indicate  rejection 
at  low  significance  levels  for  all  fifteen  series.   The  results  imply 
that  there  is  some  persistence  in  volatility.   With  the  parameter 
estimates,  we  now  turn  to  the  pricing  of  options. 
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III.   Empirical  Tests  of  the  Model 

In  this  section  we  test  the  random  variance  model  to  see  how  well 
it  prices  call  options.   We  compare  the  random  variance  model  with  the 
Black-Scholes  model  and  examine  some  of  the  pricing  biases  which  have 
been  documented  in  the  option  pricing  literature.   Although  the  Black- 
Scholes  model  is  based  on  volatility  remaining  constant,  we  apply  the 
model  with  an  ISD  that  changes  each  day.   This  procedure  has  been  used 
frequently  in  the  option  pricing  literature.   Our  data  consist  of 
transaction  prices  for  call  options  traded  on  the  CBOE ;  the  data  are 
available  in  the  Berkeley  Options  database.   We  use  prices  for  the 
year  1983.   The  Berkeley  database  contains  synchronous  stock  prices, 
but  it  does  not  contain  interest  rates.   Interest  rates  on  Treasury 
bills  were  collected  from  the  Wall  Street  Journal  each  day  for  the  ten 
maturities  closest  to  two  weeks,  one  month,  two  months,  ...,  eight 
months,  and  nine  months.   For  each  option  on  each  day  we  interpolate 
to  compute  the  relevant  interest  rate. 

For  both  models,  we  need  to  compute  volatility  estimates  each  day. 
We  use  transaction  prices  for  options  with  strike  prices  that  are 
within  $10  of  the  current  price  to  compute  ISD's.   To  avoid  the  prob- 
lem of  estimating  ISD's  and  pricing  options  simultaneously,  we  use 
prices  from  11:00  a.m.  to  12:00  noon  to  compute  ISD's  and  then  price 
options  that  trade  between  1:00  and  2:30  p.m.   Our  procedure  for  com- 
puting the  ISD  is  one  that  is  commonly  used  in  the  literature:   we  set 
the  ISD  to  minimize  the  sum  of  squared  errors  between  the  model  prices 
and  actual  option  prices.   Estimating  the  ISD  for  the  Black-Scholes 
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model  is  a  well-known  procedure,  but  this  is  not  the  case  for  the  ran- 
dom variance  model.   We  use  a  Monte  Carlo  simulation  of  the  model  in 
equation  (4)  to  compute  prices  for  the  random  variance  model  and  the 
computations  proceed  as  follows: 

H(P  ,at,t,T)  -^  Z      [PtN(du)  -  Xe"r(T_t)N(d2i)]       (5) 
i=l 

ln(Pt/X)  +  r(T-t)  +yV. 

where    d.  .  ■  

11  AT 

i 

2i    li     i 

T   2 

and  V.  is  the  value  of  V  =  (  a,    Nds  for  the  i'th  simulation  and  N  is 
i  ;    (s) 

2 
the  number  of  simulations  used,   a,  x  is  simulated  for  daily  incre- 

(s) 

ments.  a     is  the  current  value  of  the  ISD  and  it  is  the  starting 

2 
point  for  our  simulation  of  a,    N  and  V.   Our  search  procedure  for  a 

(s)  t 

uses  both  first  and  second  partial  derivatives  of  H  with  respect  to 

a     so  that  we  reduce  the  number  of  iterations  necessary  for  conver- 

t  ' 

gence.   We  find  that  1,000  simulations  with  the  antithetical  variate 
method  are  sufficient  to  produce  accurate  model  prices:   the  largest 
95%  confidence  intervals  are  +lV2  cents.   Because  the  estimation  of  the 
ISD  in  the  random  variance  model  is  very  slow  and  expensive  on  a  typi- 
cal mainframe  computer,  we  have  used  a  Cray  X-MP  supercomputer.   After 
the  ISD  has  been  estimated,  simulations  to  price  options  can  be  done 
on  a  standard  computer.   We  have  performed  all  of  the  calculations  in 
this  paper  on  the  Cray  supercomputer. 
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In  the  random  variance  model,  we  have  one  additional  parameter  to 
estimate:   the  risk  premium  associated  with  volatility.   This  parameter 
cannot  be  estimated  directly  from  price  changes  on  the  underlying 
stock.   We  instead  infer  an  estimate  from  option  prices  by  using  a 
procedure  that  is  similar  to  the  estimation  of  the  ISD's.   We  use  the 
first  two  weeks  of  prices  in  January  1983  to  estimate  an  implied  value 
of  the  risk  premium  for  each  stock.   If  we  treat  the  volatility  risk 
premium,  X ,  as  a  constant,  we  have  the  following  model  for  the  risk- 
adjusted  volatility  process: 

dZna  =  [3(a-£na)-X*]dt  +  ydz  . 

The  corresponding  discrete  approximation  is 

£na  =  (a-X)  +  plno     ,  +  e  . 

The  risk  premium  is  effectively  an  adjustment  to  the  mean  of  our 
estimated  volatility  process.   Using  the  first  two  weeks  of  prices  as 
a  preliminary  estimation  period,  we  set  X,  as  well  as  the  daily  ISD's, 
to  minimize  the  squared  error  between  actual  prices  and  model  prices. 
We  then  use  that  X  estimate  to  compute  ISD's  and  price  options  for 
the  remainder  of  the  year.   The  X  adjustments  for  each  stock  are  pre- 
sented in  Table  II.   The  X  adjustments  that  we  estimate  here  can  be 
viewed  alternatively  as  adjustments  on  the  mean  reverting  level  due 
to  error  in  our  estimation  of  the  parameters  in  the  volatility  pro- 
cess.  We  note  that  in  Table  II,  we  have  an  upward  adjustment  from  a 
to  (a-X)  for  nine  of  the  thirteen  stocks.   An  upward  adjustment  cor- 
responds to  a  negative  risk  premium  for  volatility.   For  Exxon,  we  are 
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unable  to  compute  the  ISD  in  the  random  variance  model  for  subsequent 
days  if  we  use  the  X  adjustment.   We  use  the  original  a  estimate  for 
Exxon,  but  we  use  (a-X)  for  all  of  the  other  stocks. 

Before  we  turn  to  a  discussion  of  the  results,  we  need  to  describe 
our  method  of  adjusting  for  cash  dividends.   The  Dow  Jones  companies 
all  pay  cash  dividends  and  we  use  the  standard  method  for  handling 
dividends.   We  calculate  two  prices  for  each  option.   First,  we  price 
the  option  as  if  it  were  exercised  at  maturity  and  we  subtract  the 
present  value  of  cash  dividends  from  the  current  stock,  price.   Second, 
we  price  the  option  as  if  it  were  exercised  one  day  before  the  next 
ex-dividend  date,  and  we  take  the  higher  of  the  two  prices.   Geske  and 
Roll  (1984)  and  Whaley  (1982)  have  argued  that  this  method  introduces 
a  bias  in  the  pricing  of  options.   Whaley  has  proposed  a  method  for 
handling  options  on  stocks  that  have  one  cash  dividend  before  the 
expiration  of  the  option,  but  many  of  the  options  in  our  sample  have 
more  than  one  cash  dividend  coming  before  expiration.   To  assess  any 
differences  that  might  arise  due  to  dividends,  we  have  included  in  our 
sample  a  stock,  National  Semiconductor,  that  does  not  pay  cash  divi- 
dends.  We  also  note  that  options  on  National  Semiconductor  trade  very 
actively. 

Given  our  parameter  estimates  and  the  daily  values  of  the  ISD,  we 
price  options  for  the  thirteen  stocks  from  the  middle  of  January  to 
the  end  of  December  1983.  We  first  examine  the  ability  of  each  model 
to  fit  actual  option  prices.  The  results  are  contained  in  Table  III, 
where  we  report  root  mean  squared  errors  and  mean  absolute  errors  be- 
tween model  prices  and  actual  prices.   Overall,  the  Black-Scholes 
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model  performs  marginally  better:   both  the  root  mean  squared  error 
and  the  mean  absolute  error  are  lower  by  2  cents  on  the  98,699  options 
in  our  sample.   For  nine  of  the  thirteen  stocks,  the  Black-Scholes 
model  outperforms  the  random  variance  model.   For  four  of  the  stocks, 
including  National  Semiconductor  which  does  not  pay  cash  dividends, 
the  random  variance  model  performs  marginally  better. 

In  Tables  IV,  V,  and  VI,  we  present  regression  tests  for  biases 
in  both  option  pricing  models.   Table  IV  contains  the  regression  re- 
sults for  the  strike  price  bias:   we  regress  the  percentage  pricing 
error  against  the  percent  in  or  out-of -the-money.   The  percentage 
pricing  error  is 

C,   -  H.  (a.  ) 

m     it    it   it 

yit   "   H.  (a.  ) 
it   it 

and  the  strike  price  measure  is 

P.   -X.e-r«-C) 

It        1 

"it  '    x  e-r(T-t)    • 

1 

For  our  regression  tests,  we  use  all  of  the  options  for  a  given 
stock.   The  slope  coefficients  for  the  Black-Scholes  model  are  nega- 
tive for  all  13  stocks  and  11  of  these  are  statistically  significant 
at  the  5%  level.   A  negative  coefficient  implies  that  the  model  tends 
to  underprice  out -of -the-money  calls  and  overprice  in-the-money  calls. 
The  results  for  the  random  variance  model  are  mixed.   Three  of  the 
slope  coefficients  are  positive  and  two  of  these  are  significant. 
Nine  of  the  ten  negative  coefficients  are  significant.   In  all  but 
two  of  the  cases,  the  negative  t  statistics  are  smaller  in  absolute 
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value  for  the  random  variance  model.   These  results  indicate  a  con- 
sistent negative  strike  price  bias  for  the  Black-Scholes  model.   There 
is  also  evidence  of  pricing  biases  with  the  random  variance  model,  but 
we  see  both  positive  and  negative  biases.   Geske  and  Roll  (1984)  have 
attributed  some  of  this  strike  price  bias  to  the  payment  of  cash  divi- 
dends, but  we  also  find  a  strong  negative  bias  for  National  Semi- 
conductor, which  does  not  pay  cash  dividends. 

Table  V  contains  the  regression  results  for  the  time  to  maturity 
bias.   Again  the  slope  coefficients  are  all  negative  for  the  Black- 
Scholes  model,  but  two  of  these  are  not  statistically  significant. 
The  slope  coefficient  for  GM  is  very  close  to  zero  for  the  Black- 
Scholes  model.   This  negative  bias  for  the  Black-Scholes  model  implies 
a  tendency  to  underprice  short  time  to  maturity  options  and  overprice 
long  time  to  maturity  options.   Again  the  results  are  mixed  for  the 
random  variance  model:   four  of  the  slope  coefficients  are  positive 
and  nine  are  negative.   Eleven  of  these  thirteen  coefficients  are 
significantly  different  from  zero.   The  results  in  Table  VI  for  the 
volatility  bias  are  mixed  for  both  models.   Four  of  the  thirteen  slope 
coefficients  are  positive  for  the  Black-Scholes  model  and  two  of  the 
thirteen  coefficients  are  positive  for  the  random  variance  model.   Ten 
of  the  slope  coefficients  are  statistically  different  from  zero  for 
the  Black-Scholes  model  and  nine  are  statistically  different  from  zero 
for  the  random  variance  model.   A  negative  coefficient  in  Table  VI 
implies  that  the  model  underprices  options  when  volatility  is  low  and 
overprices  options  when  volatility  is  high. 
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The  three  tables,  IV,  V,  and  VI,  indicate  evidence  of  biases  in 
both  models  with  respect  to  actual  prices.   There  is  clear  evidence 
that  the  Black-Scholes  model  has  negative  biases  with  respect  to  the 
strike  price  and  time  to  maturity.   The  results  for  the  random  var- 
iance model,  by  contrast,  are  mixed.   In  Table  VII,  we  present  the 
pricing  errors  by  nine  different  cells.   Here  we  consider  both  the 
strike  price  measure  and  time  to  maturity.   Within  each  cell  we  cal- 
culate the  average  of  the  percentage  pricing  errors  and  the  correspond- 
ing t  statistics.   Six  of  the  nine  cells  indicate  a  significant  bias 
for  the  random  variance  model,  and  eight  of  the  nine  cells  indicate 
a  significant  bias  for  the  Black-Scholes  model.   For  example,  both 
models  underprice  out-of-the-money  options  with  less  than  92  days  to 
expiration.   The  negative  strike  price  bias  found  in  Table  IV  appears 
to  be  largely  the  result  of  the  large  underpricing  errors  for  out-of- 
the-money  options  with  less  than  91  days  to  maturity.   Both  models 
overprice  in-the-money  options  with  short  maturities,  but  they  under- 
price in-the-money  options  with  medium  and  long  maturities  (over  91 
days).   Even  though  the  regression  tests  indicate  a  negative  bias  for 
tirae-to-maturity ,  the  cell  analysis  indicates  that  the  random  variance 
model  underprices  options  with  long  maturities  (over  182  days).   Our 
cell  analysis  confirms  the  evidence  of  biases  found  with  the  regression 
tests,  but  a  different  pattern  emerges.   The  random  variance  model 
underprices  long  time-to-maturity  options  and  out-of-the-money  options 
with  short  maturities;  it  appears  to  overprice  only  in-the-money 
options  with  short  maturities.   The  Black-Scholes  model  underprices 
out-of-the-money  options  with  short  maturities  and  in-the-money 
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options  with  medium  and  long  maturities;  it  overprices  options  in  the 
other  cells.   The  results  presented  here  on  the  strike  price  bias  for 
the  Black-Scholes  model  in  1983  are  exactly  opposite  to  those  reported 
by  MacBeth  and  Merville  (1979).   Rubinstein  (1985)  has  noted  that 
there  have  been  reversals  in  the  strike  price  bias  for  the  Black- 
Scholes  model. 

IV.   Hedging  Volatility  Changes 

We  turn  now  to  the  hedging  of  options  with  a  random  variance 
model.   It  has  been  noted  in  the  previous  work  on  random  variance 
option  pricing  that  one  cannot  create  a  riskless  hedge  with  the  stock 
and  an  option  as  in  the  Black-Scholes  model.   The  problem  that  arises 
is  the  extra  source  of  random  variation  in  the  option  due  to  volatil- 
ity changes.   In  the  random  variance  model,  one  can  create  a  dynamic 
portfolio  that  is  risk-free  by  hedging  an  option  position  with  the 
stock  and  another  option.   Traders  routinely  apply  the  hedging  prin- 
ciple of  the  Black-Scholes  model  by  computing  a  delta,  the  partial 
derivative  of  the  option  pricing  function  with  respect  to  the  stock 
price,  and  hedging  one  option  with  another  option.   In  the  random 
variance  model,  one  must  use  two  other  options  in  order  to  hedge  both 
the  price  risk  and  the  volatility  risk  of  an  option  position.   We  test 
this  implication  of  the  random  variance  model  by  computing  the  vari- 
ability for  hedged  positions  in  our  sample  of  options. 

Following  Hull  and  White,  we  refer  to  the  hedging  of  price  and 
volatility  risk  as  delta-slgma  hedging.  The  hedging  technique  em- 
ployed by  most  traders  is  delta  hedging.   The  delta-sigma  hedge  for  a 
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call  option  is  formed  by  creating  the  following  dynamic  portfolio 
which  is  rebalanced  at  regular  time  intervals: 


Cl  +  W2C2  +  W3C3> 

&1V3  '  Vl  '(V2  '  Vl? 

where    w„  = ; ,  w„  ■  —7- : 7  , 

2   A3v2  -  A2v3  3    (A3v2  -  A2v3) 

3C  3C. 

Ai  =  IF  '     and  Vi  =  J^   ' 

We  compute  A.  and  v.  from  our  option  pricing  functions.   If  we  use  the 

stock  in  place  of  our  third  option  we  have 


Cl  +  W2C2  +  V, 


—v  V 

1  1 

where  w0  =  and  w0  -  -A-  +  A„  — .   We  use  the  second  option  to 

2    v2      3     1    2  v2 

hedge  volatility  risk  and  the  stock  to  hedge  price  risk  in  both  op- 

Al      X 

tions.   For  delta  hedging  we  use  (C   -  A  P)  or  (C   — — -  C„). 

We  use  both  the  Black-Scholes  and  the  random  variance  models  to 
compute  the  derivatives  and  hedge  ratios  for  our  sample.   We  find  that 
the  A's  are  virtually  the  same  for  both  models,  but  the  v's  and  the 
hedge  ratios  for  delta-sigraa  hedging  are  quite  different.   For  each 
stock  in  our  sample  we  hedge  an  at-the-money  option  for  one  trading 
day.   We  revise  our  hedged  portfolios  each  day.   We  find  that  only 
eight  of  the  thirteen  stocks  have  sufficient  trading  activity  in  their 
options  to  compute  gains  on  hedged  portfolios  for  every  day  in  our 
sample.   These  stocks  are  shown  in  Tables  VIII  and  IX.   The  net  gain 
for  each  zero-investment  hedged  position  is  computed  as  follows: 
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net  gain  -  (C^  -C     )  +  w2(C2t  -C     ) 


+  w,(C_   -  C,   .)  -  r_(C      +  w0C.     +  w,C,   .), 
3   3t    3,t-l     f   l,t-l    2  2,t-l    3  3,t-l 

where  r,  *-s  c^e  risk-free  rate  calculated  for  the  holding  period. 


When  we  hedge  with  stock,  we  use  the  stock  in  place  of  C  .   For  delta 

Al 
hedging  we  set  w~  -  0  and  w  -  -  — . 

We  compute  the  normalized  net  gain  which  is  the  net  gain  above 

divided  by  C     ,  the  price  of  the  option  being  hedged.   We  have  per- 
J      l,t-l 

formed  two  sets  of  calculations.   First,  we  compute  the  standard 
deviation  for  the  normalized  net  gain  over  intervals  of  one  trading 
day.   Some  of  the  variation  in  the  transaction  prices  is  the  result  of 
variation  within  the  bid-ask  spread.   To  assess  the  impact  of  this 
extra  noise,  we  recalculate  the  standard  deviation  for  hedged  port- 
folios over  intervals  of  two  trading  days.   The  results  for  one-day 
intervals  are  contained  in  Table  VIII  and  the  results  for  two-day 
intervals  are  contained  in  Table  IX.   We  note  that  only  seven  of  the 
stocks  have  been  used  in  the  hedging  analysis:   the  other  stocks  did 
not  have  all  of  the  transaction  prices  between  1:00  and  2:30  p.m. 
needed  for  computing  gains  on  hedge'd  portfolios.   In  all  cases,  delta 
hedging  represents  a  substantial  decrease  in  variability  versus  no 
hedge.   If  we  are  hedging  with  the  underlying  stock  over  one-day 
intervals,  the  results  are  rather  poor  for  delta-sigma  hedging.   If  we 
exclude  GM,  the  variability  is  higher  for  delta-sigma  hedging  with  the 
random  variance  model.   For  four  of  the  eight  stocks,  variability  is 
lower  for  delta-sigma  hedging  with  the  Black -Scholes  model.   If  we  are 
hedging  with  options  only  over  one-day  intervals,  the  results  for 
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delta-sigma  hedging  are  much  better:   variability  with  delta-sigraa 
hedging  is  lower  for  both  models.   The  two  exceptions  are  Exxon  with 
the  random  variance  model  and  AT&T  with  the  Black-Scholes  model.   The 
results  are  quite  different  in  Table  IX  for  two-day  intervals:   the 
variability  is  lower  with  delta-sigma  hedging  for  most  stocks.   In 
eight  of  the  fourteen  comparisons  for  the  random  variance  model,  the 
variability  of  delta-sigma  hedging  is  lower  than  that  of  delta  hedging. 
In  ten  of  the  fourteen  comparisons  for  the  Black-Scholes  model,  the 
variability  of  delta-sigma  hedging  is  lower.   In  some  cases,  the 
reduction  in  the  standard  deviation  is  as  much  as  20  to  30  percent. 
The  results  in  Table  IX  provide  support  for  delta-sigma  hedging  as  a 
means  of  reducing  risk  with  stock  options.   Hull  and  White  (1987b)  and 
Chesney  and  Scott  (1988)  have  found  that  delta-sigraa  hedging  reduces 
substantially  the  variability  of  hedged  positions  in  foreign  currency 
options. 

V.   Summary  and  Conclusions 

In  this  paper  we  have  presented  an  empirical  analysis  of  the  ran- 
dom variance  option  pricing  model  by  using  actual  transaction  prices 
on  call  options  from  the  CBOE  for  the  year  1983.   By  using  transaction 
prices  and  synchronous  stock  prices,  we  have  avoided  some  of  the  pit- 
falls of  earlier  empirical  studies.   To  implement  the  random  variance 
model,  we  use  parameter  values  for  the  volatility  process  which  are 
estimated  from  stock  return  data.   We  use  both  the  random  variance 
model  and  the  Black-Scholes  model  to  compute  ISD's  and  to  price  a  very 
large  number  of  options.   Both  models  are  able  to  explain  the  observed 
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prices  and  the  Black-Scholes  model  performs  marginally  better  than  the 
random  variance  model.   The  difference  in  pricing  errors  between  the 
models  averages  2  cents. 

We  find  that  there  are  consistent  biases  in  pricing  by  the  Black- 
Scholes  model  with  respect  to  the  strike  price  and  time  to  maturity. 
We  also  find  biases  with  the  random  variance  model,  but  the  biases  are 
mixed.   Some  are  positive  and  some  are  negative.   We  include,  in  our 
sample,  options  on  a  stock  that  does  not  pay  cash  dividends  and  we 
find  that  the  results  are  similar  to  those  for  options  on  stocks  that 
do  pay  cash  dividends.   As  in  Rubinstein's  (1985)  study  of  options 
from  1976  to  1978,  we  find  that  cash  dividends  are  not  an  explanation 
for  these  observed  biases.   In  the  last  section  of  the  paper  we 
examine  the  benefits  of  hedging  both  price  and  volatility  risk  and  we 
find  that  one  can  reduce  the  variability  of  an  option  position  by 
hedging  volatility  risk. 
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FOOTNOTES 

1.  To  get  Che  fundamental  P.D.E.,  we  use  Theorem  3  of  Cox,  Ingersoll, 
and  Ross  (1985a)  and  the  observation  in  their  companion  paper, 
CIR  (1985b),  that  wealth  drops  out  of  the  model  if  we  have  con- 
stant relative  risk  aversion  utility.   In  our  model,  we  have  only 
two  random  state  variables,  the  asset  price  and  volatility.   Even 
if  we  have  more  state. variables ,  for  a  standard  call  option  written 
on  a  stock  we  can  reduce  the  option  pricing  function  to  a  function 
of  time,  price,  and  volatility.   The  solution  to  the  fundamental 
P.D.E.  follows  from  lemma  4  of  CIR. 

2.  Wiggins  has  computed  estimates  for  the  correlation  between  stock 
price  changes  and  volatility  and  he  also  finds  that  most  of  his 
estimates  are  not  statistically  different  from  zero. 

3.  A  discussion  of  approximations  over  short  time  intervals  can  be 
found  in  Arnold  (1974,  pp.  39-41)  and  Karlin  and  Taylor  (1981, 
pp.  340-58).   This  approximation  for  A£nP  can  be  obtained  by 
replacing  the  differentials  in  dilnP  with  differences.   We  also 
note  that  Monte  Carlo  simulation  of  a  stochastic  volatility 
option  pricing  model  will  require  simulations  over  small  discrete 
time  intervals. 

4.  Note  that  we  have  used  the  following  result:   E(e    ta|u    |)  =  0 

even  though  E(e   u   )  =  8a   .   We  use  the  result  again  later  in 
t-1  t-l      e 

our  analysis.   The  tedious  mathematical  results  used  in  this 


-25- 

section  are  contained  in  an  appendix  which  can  be  obtained  from 
the  author. 

5.  For  a  reference,  see  Nerlove,  Grether,  and  Carvalho  (1979),  pp. 
97-102). 

6.  See  Dhrymes  (1970,  p.  113). 

7.  Each  one  of  the  sample  moments  is  a  function  of  the  data:   m.  = 
—  Z  f  (x  )•   Consider  a  vector  of  these  time  series  with  the  ith 

T  t-i  1  c 

element,  (f.(x  )  -  m.),  and  the  corresponding  spectral  density 
matrix.   The  variance  matrix  for  the  sample  moments  is  equal  to  2tt 
times  the  spectral  density  matrix  evaluated  at  the  zero  frequency 
divided  by  the  sample  size,  T.   If  x  were  independent  over  time 
and  all  of  our  functions  of  the  data,  f.(x  ),  were  serially  uncor- 
rected, then  the  variance  matrix  for  the  sample  moments  would 
collapse  to  the  familiar  result  which  is  the  variance  matrix  for 
the  time  series  divided  by  T.   In  our  model,  the  vector  of  time 
series  is  serially  correlated  and  we  allow  for  this  by  using  the 
spectral  estimator.   To  estimate  the  spectral  density  matrix  at  the 
zero  frequency,  we  use  the  smoothed  periodogram  estimator  and  the 
procedure  described  in  Nerlove,  Grether,  and  Carvalho  (1979,  pp. 
67-68).   We  prewhiten  the  series  and  estimate  the  spectrum  of  the 
prewhitened  series  with  a  flat  window.   The  estimates  are  then  re- 
colored  by  the  appropriate  transfer  functions. 


8.   For  x  series  that  we  use  is  also  the  residual  from  a  dummy  var- 
iable regression. 
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TABLE  II 

Risk  Premia  Adjustments  Estimated 
From  Option  Prices 


ina     =  (a-X*)  +  plna         +   e 


Implied 

Mean  Reverting 

Adjusted 

Mean  Reverting 

Stock 

a 

a   Level 

(a-X) 
-.02586 

a  Level 

ALCOA 

-.02632 

.01511 

.01627 

AT&T 

-.01997 

.00733 

-.01633 

.01795 

Exxon 

-.14653 

.01093 

-.13296* 

.01660 

GE 

-.07818 

.01186 

-.07694 

.01272 

GM 

-.01123 

.01300 

-.01225 

.00876 

IBM 

-.01670 

.01117 

-.01798 

.00791 

Inter. 

Paper 

-.04437 

.01404 

-.04484 

.01342 

Kodak 

-.05717 

.01459 

-.05646 

.01537 

Merck 

-.22053 

.01274 

-.21517 

.01417 

MMM 

-.04000 

.01273 

-.03855 

.01492 

N.  Semi. 

-.02400 

.02905 

-.02712 

.01833 

Sears 

-.03672 

.01381 

-.03384 

.01933 

United 

Tech. 

-.02361 

.01502 

-.02334 

.01577 

NOTE:   a  level  is  the  daily  standard  deviation. 

a     (a-X) 

The  mean  reverting  level  is  ~ or  — 

a  1-p     1-p 


*With  the  X  adjustment  on  Exxon,  we  were  unable  to  estimate  the  ISD 
for  subsequent  days.   For  this  reason,  we  have  not  used  the  X  adjust- 
ment on  Exxon. 


TABLE 

III 

Model  Comparison 

Random 

Variance 

Black- 

Scholes 

Model 

Model 

Root  Mean    Mean 

Root  Mean 

Mean 

Number 

Squared 

Absolute 

Squared 

Absolute 

Stock 

of  Options 

Error 

Error 

Error 

Error 

ALCOA 

2,385 

$.2419 

$.1849 

$.2433 

$.1850 

AT&T 

8,947 

.2954 

.2131 

.2706 

.1957 

Exxon 

5,338   • 

.3904 

.1326 

.3791 

.1091 

GE 

10,372 

.3260 

.2290 

.2918 

.1984 

GM 

12,195 

.2767 

.2002 

.2280 

.1611 

IBM 

21,135 

.4132 

.2497 

.3876 

.2149 

Inter.  Paper 

1,639 

.2708 

.1986 

.2716 

.1995 

Kodak 

11,262 

.3061 

.1654 

.2990 

.1586 

Merck 

2,071 

.5017 

.3501 

.3830 

.2634 

MMM 

2,919 

.3005 

.2105 

.2409 

.1631 

N.  Semi. 

9,781 

.2685 

.1839 

.3092 

.2147 

Sears 

8,481 

.1906 

.1247 

.1851 

.1199 

United  Tech. 

2,174 

.2632 

.1957 

.2679 

.1972 

All  Options 

98,699 

.3262 

.2022 

.3060 

.1825 

TABLE  IV 

Regression  Tests 
Strike  Price  Bias 


y   =  a  +  bm   +  e 
3  it         it    it 


Stock 

R2 

(Number  of  Options) 

Model 
RV 

a 

t(a=0) 
5.21 

b 

t(b=0) 
-6.96 

ALCOA 

.0184 

-.2085 

.0199 

(2,385) 

BS 

.0296 

7.66 

-.2824 

-8.62 

.0302 

AT&T 

RV 

.6901 

4.90 

-11.7549 

-6.14 

.0042 

(8,947) 

BS 

.1570 

5.52 

-2.7888 

-7.20 

.0058 

Exxon 

RV 

.3116 

3.50 

-2.5690 

-3.42 

.0022 

(5,337) 

BS 

.3118 

4.83 

-2.7090 

-4.98 

.0046 

General  Electric 

RV 

.0686 

6.47 

-.7253 

-7.45 

.0053 

(10,372) 

BS 

.1865 

7.07 

-2.0151 

-8.33 

.0066 

General  Motors 

RV 

-.0294 

-8.53 

.1861 

5.48 

.0025 

(12,195) 

BS 

.0001 

.01 

-.0583 

-1.03 

.0001 

IBM 

RV 

-.0221 

-43.51 

.1157 

30.21 

.0414 

(21,135) 

BS 

.0028 

5.11 

-.0229 

-5.48 

.0014 

International  Paper 

RV 

-.0110 

-2.90 

.0770 

1.92 

.0022 

(1,638) 

BS 

.0063 

1.20 

-.1034 

-1.87 

.0021 

Kodak 

RV 

.1783 

1.98 

-6.3181 

-5.57 

.0027 

(11,262) 

BS 

.2700 

1.99 

-16.5616 

-9.66 

.0082 

3M 

RV 

.0805 

11.05 

-.7289 

-8.21 

.0226 

(2,919) 

BS 

.0293 

7.98 

-.2857 

-6.37 

.0137 

Merck 

RV 

1.7894 

2.20 

-18.0846 

-1.75 

.0015 

(2,071) 

BS 

.0834 

8.68 

-1.3428 

-11.00 

.0052 

National  Semiconductor 

RV 

.0015 

1.33 

-.0194 

-3.80 

.0015 

(9,781) 

BS 

.0381 

14.93 

-.1232 

-10.63 

.0114 

Sears 

RV 

.1200 

4.39 

-.7817 

-4.23 

.0021 

(8,481) 

BS 

.1000 

5.11 

-.6820 

-5.16 

.0031 

United  Technology 

RV 

.0664 

11.04 

-.4573 

-7.59 

.0258 

(2,174) 

BS 

.0554 

10.85 

-.3867 

-7.57 

.0257 

TABLE  V 

Regression  Tests 
Time-to-Maturity  Bias 


v .   =  a  +  dt.   +  e. 
'it         it    it 


Stock 

R2 

(Number  < 

af  Options) 

Model 
RV 

a 

t(a-0) 
3.46 

b 

t(b=0) 
-1.95 

ALCOA 

.0232 

-.0001 

.0016 

(2,385) 

BS 

.0539 

7.35 

-.0003 

-5.25 

.0114 

AT&T 

RV 

1.6535 

6.01 

-.0097 

-4.94 

.0027 

(8,947) 

BS 

.3737 

6.71 

-.0022 

-5.55 

.0034 

Exxon 

RV 

.5659 

3.34 

-.0027 

-2.46 

.0011 

(5,337) 

BS 

.6252 

5.11 

-.0032 

-4.01 

.0030 

General  Electric 

RV 

.0630 

3.57 

-.0002 

-1.28 

.0002 

(10,372) 

BS 

.2950 

6.73 

-.0021 

-4.89 

.0023 

General  Motors 

RV 

-.0692 

-12.37 

.0005 

10.23 

.0085 

(12,195) 

BS 

.0003 

.04 

-.0000 

-.33 

.0000 

IBM 

RV 

-.0482 

-63.35 

.0003 

53.35 

.1187 

(21,135) 

BS 

.0124 

14.35 

-.0001 

-15.41 

.0111 

International  Paper 

RV 

-.0265 

-3.87 

.0002 

3.07 

.0057 

(1,638) 

BS 

.0144 

1.52 

-.0001 

-1.40 

.0012 

Kodak 

RV 

.8347 

5.05 

-.0060 

-.407 

.0015 

(11,262) 

BS 

1.4377 

5.74 

-.0098 

-4.42 

.0017 

3M 

RV 

.1895 

15.93 

-.0016 

-13.07 

.0553 

(2,919) 

BS 

.0723 

11.95 

-.0006 

-10.13 

.0339 

Merck 

RV 

4.0517 

2.96 

-.0333 

-2.37 

.0027 

(2,071) 

BS 

.1444 

8.77 

-.0011 

-6.73 

.0214 

National 

Semiconductor 

RV 

-.0010 

-.56 

.0000 

.61 

.0000 

(9,781) 

BS 

.0909 

22.31 

-.0007 

-18.87 

.0351 

Sears 

RV 

.2694 

5.58 

-.0018 

-4.56 

.0024 

(8,481) 

BS 

.2088 

6.04 

-.0013 

-4.82 

.0027 

United  Technology 

RV 

.1247 

13.38 

-.0009 

-10.07 

.0445 

(2,174) 

BS 

.0956 

12.03 

-.0007 

-8.57 

.0327 

TABLE  VI 

Regression  Tests 
Volatility  Bias 


y.   =  a  +  bo.      +  e. 
'it         it    it 


Stock 

R2 

(Number  < 

jf  Options) 

Model 
RV 

a 

t(a=0) 
10.23 

b 

t(b=0) 
-9.86 

ALCOA 

.2671 

-12.4211 

.0391 

(2,385) 

BS 

.2668 

7.20 

-12.2347 

-6.67 

.0183 

AT&T 

RV 

6.3110 

6.11 

-485.4491 

-5.70 

.0036 

(8,947) 

BS 

1.8256 

6.99 

-125.9132 

-6.63 

.0049 

Exxon 

RV 

.7057 

2.47 

-33.6802 

-1.84 

.0006 

(5,337) 

BS 

.6175 

1.10 

-32.1893 

-.75 

.0001 

General  i 

Electric 

RV 

.1486 

2.58 

-5.6129 

-1.83 

.0003 

(10,372) 

BS 

.0745 

.33 

11.4779 

.86 

.0001 

General  Motors 

RV 

.1161 

3.87 

-7.3080 

-4.62 

.0018 

(12,195) 

BS 

.1535 

2.72 

-8.7553 

-2.77 

.0006 

IBM 

RV 

-.0330 

-9.37 

1.0365 

5.32 

.0013 

(21,135) 

BS 

-.0579 

-13.18 

3.5796 

13.56 

.0086 

International  Paper 

RV 

.0177 

7.18 

-10.8689 

-7.61 

.0341 

(1,638) 

BS 

.1530 

3.04 

-8.9677 

-2.99 

.0054 

Kodak 

RV 

.2611 

.50 

.3772 

.01 

.0000 

(11,262) 

BS 

-2.3489 

-1.82 

173.2574 

2.23 

.0004 

3M 

RV 

.3439 

7.95 

-23.1840 

-6.60 

.0147 

(2,919) 

BS 

.1525 

4.71 

-9.7688 

-4.05 

.0056 

Merck 

RV 

6.8834 

2.95 

-534.1532 

-2.50 

.0030 

(2,071) 

BS 

.3968 

3.88 

-25.0538 

-3.37 

.0055 

National 

Semiconductor 

RV 

.0969 

10.79 

-3.0073 

-10.88 

.0119 

(9,781) 

BS 

.0994 

4.05 

-2.4319 

-2.92 

.0009 

Sears 

RV 

.2379 

1.14 

-7.7331 

-.74 

.0001 

(8,481) 

BS 

.0258 

.12 

2.052 

.19 

.0000 

United  Technology 

RV 

.3123 

7.39 

-17.0729 

-6.29 

.0178 

(2,174) 

BS 

.3374 

7.87 

-18.7336 

-6.96 

.0218 

TABLE  VII 

Cell  Analysis  of  Option  Pricing  Biases 
Average  Percentage  Pricing  Errors  and  t  Statistics 


Strike  Price 

out -of -the -money 


at-the-money 


in-the-money 


Time  to  Maturity 


Short 


10,404 


U 


RV   1.1559  (5.20) 


23,642 

RV  -.0014  (-1.00) 


20,926 


RV  -.0062  (-22.10) 


BS    .9859  (6.43) 


BS   -.0021  (-2.40) 


BS   -.0010   (-3.56) 


Medium 


6,465 

RV   -.0007  (-.35) 

0 
BS   -.0035  (-1.96) 


13,750 

RV   -.0001  (-.10) 

0 
BS   -.0075  (-13.38) 


10,882 


U 


RV    .0031   (3.77) 

U 
BS    .0022   (2.72) 


Long 


1,619 


RV   .0450  (10.71) 


5,701 


RV    .0087  (6.73) 


5,308 


RV    .0150  (18.43) 


BS   -.0022  (-.71) 


BS   -.0069  (-6.78) 


BS    .0057  (7.94) 


NOTES:   For  at-the-money  options 


P   -  Xe 
t 


-r(T-t) 


Xe 


-r(T-t) 


<  .05 


Time  to  Maturity:   Short,  _<  91  days 

Medium,  92-182  days 
Long,  _>  183  days 

RV  =  random  variance  model 

BS  =  Black-Scholes  model 

U  =  underprices 

0  =  overprices 


TABLE  VIII 

Standard  Deviation  of  Normalized  Net  Gain 
One-Day  Intervals 


Random  Variance  Model 


Hedge  with  Stock 


Hedge  with  Options  Only 


Stock 

No  Hedge 

Delta  Hedge 

Delta-Sigma 

Delta  Hedge 

Delta-Sigma 

AT&T 

.1376 

.0777 

.0875 

.0989 

.0984 

Exxon 

.1642 

.1153 

.1365 

.1026 

.1221 

GM 

.1544 

.0538 

.0509 

.0579 

.0533 

IBM 

.1113 

.0313 

.0340 

.0382 

.0321 

Kodak 

.1445 

.0550 

.06173 

.0629 

.0622 

N.  Semi. 

.2088 

.1013 

.1392 

.1312 

.1246 

Sears 

.1595 

.0577 

.0969 

.0833 

.0778 

Black -Scholes  Model: 


Hedge  with  Stock 


Hedge  with  Options  Only 


Stock 

No  Hedge 

Delta  Hedge 

Delta-Sigma 

Delta  Hedge 

Delta-Sigma 

AT&T 

.1376 

.0802 

.0841 

.0990 

.1040 

Exxon 

.1642 

.1154 

.1033 

.1144 

.1090 

GM 

.1544 

.0537 

.0483 

.0583 

.0549 

IBM 

.1113 

.0312 

.0306 

.0383 

.0335 

Kodak 

.1445 

.0551 

.0487 

.0629 

.0646 

N.  Semi. 

.2088 

.1011 

.1273 

.1309 

.1271 

Sears 

.1595 

.0588 

.0755 

.0837 

.0732 

TABLE  IX 

Standard  Deviation  of  Normalized  Net  Gain 
Two-Day  Intervals 


Random  Variance  Model 


Hedge  with  Stock 


Hedge  with  Options  Only 


Stock 

No  Hedge 

Delta  Hedge 

Delta-Sigma 

Delta  Hedge 

Delta-Sigma 

AT&T 

.1998 

.0896 

.0788 

.1006 

.1062 

Exxon 

.2405 

.1139 

.1228 

.0999 

.1230 

GM 

.2211 

.0753 

.0627 

.0658 

.0591 

IBM 

.1602 

.0419 

.0373 

.0427 

.0357 

Kodak 

.1847 

.0722 

.0680 

.0745 

.0667 

N.  Semi. 

.3097 

.1391 

.1913 

.1721 

.1348 

Sears 

.2231 

.0750 

.1150 

.0877 

.0997 

Black-Scholes  Model 


Hedge  with  Stock 


Hedge  with  Options  Only 


Stock 

No  Hedge 

Delta  Hedge 

Delta-Sigma 

Delta  Hedge 

Delta-Sigma 

AT&T 

.1998 

.0930 

.0774 

.0997 

.1045 

Exxon 

.2405 

.1146 

.0782 

.1143 

.0921 

GM 

.2211 

.0749 

.0597 

.0654 

.0551 

IBM 

.1602 

.0413 

.0343 

.0426 

.0388 

Kodak 

.1847 

.0730 

.0572 

.0746 

.0757 

N.  Semi. 

.3097 

.1385 

.1775 

.1724 

.1697 

Sears 

.2231 

.0777 

.0925 

.0887 

.0848 
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